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The ground-state properties of spin-polarized tritium TJ at zero temperature are obtained by 
means of diffusion Monte Carlo calculations. Using an accurate ab initio TJ-TJ interatomic potential 
we have studied its liquid phase, from the spinodal point until densities above its freezing point. 
The equilibrium density of the liquid is significantly higher and the equilibrium energy of —3.664(6) 
K significantly lower than in previous approximate descriptions. The solid phase has also been 
studied for three lattices up to high pressures, and we find that hep lattice is slightly preferred. The 
liquid-solid phase transition has been determined using the double-tangent Maxwell construction; 
at zero temperature, bulk tritium freezes at a pressure of P — 9(1) bar. 



PACS numbers: 67.63.Gh,67.80.fh 

I. INTRODUCTION 

Spin-polarized hydrogen (HJ,), and its isotopes spin- 
polarized deuterium (DJ) and tritium (TJ), are fully 
quantum systems with a wealth of interesting properties. 
The low mass of HJ. and the very weak attractive part of 
its interaction potential makes possible that it remains in 
a stable gas phase even in the limit of zero temperature. 
Recent microscopic calculations at zero temperature have 
shown that the gas phase persists for pressures up to 170 
bari The gaseous state of bulk HJ. was proposed in 1976 
by Stwalley and Nosanow 2 . as the best reliable option to 
achieve a Bose-Einstein condensate (BEC) state without 
the strong depletion of the condensate that correlations 
induce in superfluid 4 He. After 22 years of intense and 
continued experimental work by many groups, Fried et 
al& managed to achieve a BEC of HJ. Problems that 
needed to be overcome included recombination on the 
walls, by working with a wall-free confinement, and low 
evaporation rates by using spin resonance. In the mean- 
time, other BEC states were first achieved in 1995 work- 
ing with alkali gases such as Rb, Na and Li A 

In contrast to alkali gases, the hydrogen-hydrogen in- 
teratomic interaction is very well known and leads to a 
stable gas phase. Spin-polarized hydrogen atoms interact 
via the triplet potential b 3 E+ determined in an essen- 
tially exact way by Kolos and Wolniewicz^ and recently 
extended to larger interparticle distances by Jamieson et 
al.& The s-wave scattering length a is appreciably smaller 
than the typical values for alkalis, which retards evapora- 
tive cooling and produces a higher transition temperature 
(50 v>K). 

Essentially the same interaction can be applied for 
heavier isotopes spin-polarized deuterium DJ and tritium 
TJ. DJ atoms obey Fermi statistics, and so the zero- 
pressure state of bulk D J depends on the number of occu- 
pied nuclear spin states. In the limit of zero pressure and 
zero temperature, previous theoretical studies£&£ have 
shown that (DJi) with only one occupied nuclear spin 



state is a gas, while bulk DJ with two (DJ.2) and three 
(DJ3) equally occupied nuclear spin states remains liq- 
uid. Spin-polarized tritium, which obeys Bose statistics, 
is expected to be liqui d 10 ' 11 ^ 2 due to its larger mass. In 
fact, Stwaley and Nosanow^ suggested it should behave 
very much like liquid 4 He and therefore constitute a sec- 
ond example of bosonic superfluid. Recently, microscopic 
properties of TJ clusters have been studied by Blume et 
al. 13 using the diffusion Monte Carlo (DMC) method. In 
addition, in Ref. [13| spin-polarized tritium is suggested 
as a new BEC gas in optical dipol trap. It has the same 
advantage of a nearly exact knowledge of the interatomic 
potential as spin-polarized hydrogen but, unlike HJ, it 
has a very broad Feshbach resonance that can be used to 
control the condensate in a trap. 

In the present work, we present a DMC study of the 
liquid and solid phases of spin-polarized tritium. For 
bosonic many-body systems at zero temperature DMC 
methods lead to exact estimations of the ground-state 
energy and related properties within statistical errors. 
Using the ab initio TJ-TJ interatomic potential within 
the DMC method, we report accurate microscopic results 
for energetic and structural properties of the bulk system. 
Relevant results of this work include the determination 
of the equilibrium density and energy per particle, the 
spinodal point and the liquid-solid phase transition. 

In Sec. II, we briefly describe the DMC method and 
the trial wave functions used for importance sampling 
in the liquid and solid phases. In Sec. Ill, the results of 
the DMC simulations are reported in several subsections. 
The first and second one are devoted to the microscopic 
results for the liquid and solid phases, respectively. In the 
last one, we study the liquid-solid phase transition point 
and report results on the freezing and melting densities. 
Finally, Sec. IV comprises a summary of our results and 
the main conclusions. 
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II. METHOD 

The starting point of the DMC method is the 
Schrodinger equation written in imaginary time, 

d^(R t) 

-h^^- = (H-E r )^(R,t) , (1) 
with an TV-particle Hamiltonian 

h 2 N N 

#=-<L£^ 2 +£^). (2) 

i— 1 i<j 

In Eq. (p}, E r is a constant acting as a reference energy 
and R = (ri, . . . , rjy) is a walker in Monte Carlo termi- 
nology. In order to reduce the variance to a manageable 
level it is a common practice to use importance sampling 
by introducing a trial wave function ip(R). Then, the 
Schrodinger equation is rewritten for the wave function 
&(R,t) = 1 i'(R } t)^/j(R) and solved in a stochastic form. 
In the limit t — > oo only the lowest energy eigenfunction, 
not orthogonal to ip(R), survives and then the sampling 
of the ground state is effectively achieved. Apart from 
statistical uncertainties, the energy of a iV-body bosonic 
system is exactly calculated. 

The interaction between TJ, atoms is described with 
the spin-independent central triplet pair potential b 3 Y,+ . 
It has been determined in an essentially exact way by 
Kolos and Wolniewicz,— and recently extended to larger 
interparticle distances by Jamieson et al. (JDW)£ As 
in a recent DMC calculation of bulk Hjji we have used 
a cubic spline to interpolate between JDW data. This 
interaction is then smoothly connected to the long-range 
behavior of the Tj-TJ, potential as calculated by Yan et 
al.t^ The JDW potential used in the present work has 
a core diameter a — 3.67 A and a minimum of —6.49 K 
at a distance 4.14 A. A comparison between different po- 
tentials employed in the past is reported in Ref. H. We 
have also verified that the addition of mass-dependent 
adiabatic corrections (as calculated by Kolos and Rych- 
lewski-i^) to the JDW potential does not change the en- 
ergy of the bulk spin-polarized tritium. 

The trial wave function used for the simulation of the 
liquid phase is of Jastrow form, 

N 

lfcr(fl) = n /fa*)- ( 3 ) 

i<j 

The two-body correlation function f(r) is the same as in 
our previous study of spin-polarized hydrogen, 1 

f(r) = exp[-6i exp(-6 2 r)] , (4) 

where b\ and 62 are variational parameters. The same 
form was also used in the variational Monte Carlo (VMC) 
calculation of Etters et al.^ who modeled the HJ,-H| in- 
teraction with a Morse potential fitted to reproduce Ko- 
los and Wolniewicz ab initio data£ This analytic form 



((4]) corresponds to the WKB solution of the two-body 
Schrodinger equation for small interparticle distances 
when the potential is of Morse type. 

Simulations of the crystalline bec, fee and hep phases 
have been also carried out; in this case, we use a 
Nosanow- Jastrow model 

AT 

^nj(JI) =MR)H Hm) , (5) 

i 

where h{r) = exp(— ar 2 /2) is a localizing function which 
links every particle i to a fixed lattice point rj. The 
parameter a is optimized variationally. 

The trial wave function ip(R) has been optimized for 
the density range where the equation of state has been 
calculated, by using the VMC method. The liquid phase 
has been studied for densities in the interval from 0.006 
A~ 3 to 0.02 A~ 3 . Within this interval, the value of the 
parameter b\ ([¥]) that optimizes the trial wave function 
takes increasing values with the density from 110 to 180, 
while the second parameter 62 Q does not change sig- 
nificantly, assuming values from 1.28 to 1.35 A -1 . For 
the three solid lattices (bec, fee and hep), the calcula- 
tions have been carried out from 0.008 to 0.024 A -3 . As 
in the case of the liquid phase, the parameter 62 slightly 
changes, from 1.29 to 1.45 A -1 . The parameter b\ in- 
creases with density from the melting point up to the 
highest density studied here, taking values from 80 to 
148. The parameter a, which models the strength of the 
localization of particles around the lattice sites, increases 
with density. In the case of the bec phase, optimized val- 
ues of a range from 0.33 to 1.80 A~ 2 , for the fee phase 
from 0.32 to 2.47 A" 2 , and for the hep phase from 0.28 
to 2.21 A -2 . The statistical errors of the variational en- 
ergies in this optimization procedure with VMC are com- 
patible with those of the DMC results (see Tables U and 
P). 

We use the DMC method accurate to second order in 
the time step At^- which allows us to use larger At val- 
ues than in linear DMC. Both the time-step dependence 
and the mean walker population have been studied care- 
fully in order to eliminate any bias coming from them. 

Any simulation of a bulk system with a finite number 
of particles requires a size-dependence analysis in order 
to achieve results as close as possible to the thermody- 
namic limit. For the liquid phase, we have used 250 par- 
ticles in all simulations and checked at the VMC level 
that with the addition of standard tail corrections, the 
size dependence of the energy remains smaller than the 
typical size of the statistical error. On the other hand, 
in all the solid state simulations we have assumed peri- 
odic boundary conditions, with 256, 250, and 180 atoms 
for the bec, fee, and hep lattices, respectively. Due to 
the periodic order of the solid, standard tail corrections 
which assume (g(r) = 1) beyond r > L/2, where L is the 
length of the simulation box, become rather inaccurate. 
In order to better determine the energy tail corrections 
we have studied the size dependence of the energy at the 
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(A- 



E/N (K) T/N (K) P (bar) c (m/s) 



0.006 

0.0074 

0.009 

0.01 

0.0125 

0.016 



-3.427(2) 
-3.664(6) 
-3.320(7) 
-2.675(8) 
0.86(2) 
12.26(6) 



5.578(17) 
7.778(18) 
10.84(3) 
13.06(4) 
19.28(6) 
30.33(10) 



-1.43(2) 
-0.12(1) 
5.11(6) 
11.6(1) 
45.2(4) 
161(2) 



71(3) 
189(3) 
318(3) 
402(4) 
630(5) 
992(8) 



TABLE I: Results for liquid TJ at different densities p: en- 
ergy per particle (E/N), kinetic energy per particle (T/N), 
pressure (P), and speed of sound (c). Figures in parenthesis 
are the statistical errors. 



VMC level, where larger number of particles can be used. 
From the VMC results one extracts the tail corrections 
for a given number of atoms and then these are added 
to the DMC energies. With this procedure it was possi- 
ble to reproduce accurately the experimental equation of 
state of solid 4 Hei£. 



III. RESULTS 
A. Liquid phase 

Spin-polarized T in its liquid phase has been studied 
in the density range from spinodal point up to densi- 
ties above crystallization. DMC results for the total and 
kinetic energy per particle at different densities are re- 
ported in Table |TJ In order to remove any residual bias 
from the trial wave function, kinetic energies are calcu- 
lated as differences between total energies and pure es- 
timations of potential energies. The total energy is neg- 
ative approximately up to the density p = 0.012 A~ 3 . 
The potential energy per particle is negative in all the 
density regime studied, presenting a minimum value of 
around -19 K at the density p = 0.014 A" 3 

In Fig. [TJ we plot the DMC results for the equation 
of state of the liquid. We have tried different analytical 
forms to fit the DMC data. The best results have been 
obtained by using a polynomial fit of the form (e = E/N) 



e(p) = e +B 



Po 



Po 



C 



P~ Po 
Po 



(6) 



Po and eo being the equilibrium density and the energy 
per particle at this density, respectively. The equation 
of state © is shown as a solid line on top of the DMC 
data in Fig. [TJ The best set of parameters is: eo = 
-3.656(4) K, B = 6.86(7) K , C = 4.70(5) K , and 
po = 0.007466(7) A -3 , the figures in parenthesis being 
the statistical uncertainties. It is worth noticing that the 
value obtained for the equilibrium density expressed in 
units of c~ 3 , pa = 0.369 er~ 3 (er = 3.67 A) is similar to 
the one in liquid 4 Hc, p = 0.365 cr" 3 (a = 2.556 A). 
Using the equation of state ©, we have obtained the 



40 
30 

Q 20 
W 10 


-10 



-3.4 
-3.6 



0.006 0.008 



0.006 0.009 0.012 0.015 0.018 0.021 
P(A" 3 ) 



FIG. 1: Energy per particle of liquid TJ. (solid circles) as a 
function of the density p. The solid line corresponds to the 
fit to the DMC energies using Eq. ((6]). The error bars of the 
DMC energies are smaller than the size of the symbols. 
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FIG. 2: Pressure and speed of sound of liquid TJ. as a func- 
tion of the density. Left (right) scale corresponds to pressure 
(speed of sound). 



pressure from its thermodynamic definition 
P(p)=p 2 (de/dp) ; 



(7) 



and from it, the corresponding speed of sound as a func- 
tion of the density 



c 2 (p) 



1 

m 



dP 
~dp~ 



(8) 



In Table [TJ we report results for the pressure P and the 
speed of sound c for some values of the density, where 
specific DMC simulations have been carried out. The 
functions P(p) and c(p), derived respectively from Eqs. 
© and ©, are shown in Fig. d 

The spinodal density in TJ, i.e., the density where the 
speed of sound becomes zero, is po = 0.0056 A -3 = 0.277 
cr~ 3 , corresponding to a pressure of P s = —1.48(2) bar. 
For comparison, the spinodal density in liquid He is a 
bit lower, po = 0.264 er -3 , the spinodal pressure being 
larger in absolute value, P s = —9.30(15) bar. In Fig. 
[3J we plot the speed of sound c ([5]) as a function of the 
pressure for pressures approximately up to solidification. 
It can be seen that c drops very fast when approaching 
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FIG. 3: Speed of sound of liquid TJ, as a function of the 
pressure from the spinodal point up to freezing. The spinodal 
pressure P s corresponds to the point where c = 0. 
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FIG. 5: Static structure function of the liquid phase. From 
bottom to top in the height of the main peak, the results 
correspond to densities 0.006, 0.0074, 0.01, 0.016, and 0.02 
A" 3 . 
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FIG. 4: Two-body radial distribution functions of the liquid 
phase. From bottom to top in the height of the main peak, 
the results correspond to densities 0.006, 0.0074, 0.01, 0.016, 
and 0.02 A" 3 . 
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FIG. 6: Condensate fraction of spin-polarized T in the liq- 
uid phase. The line corresponds to a fit to the DMC data 
using Eq.©. The error bars are smaller than the size of the 
symbols. 



the spinodal point. Near the spinodal point it is expected 
that c has the form oc (P-Ps) 1 ^, where v is the critical 
exponent. 

Apart from the ground-state energy, DMC simula- 
tions enables us to make exact estimations of other rel- 
evant magnitudes such as the two-body radial distribu- 
tion function g(r) and its Fourier transform, the static 
structure function S(k). The use of pure estimators 18 
eliminates the bias coming from the trial wave function 
and allows us to arrive to exact results for both functions. 
The evolution of g{r) with density for liquid TJ. is shown 
in Fig. [J] It is very similar to the well-known results for 
liquid 4 He. When p increases, g{r) gains structure, with 
the main peak shifting to shorter distances and increasing 
its strength. 

In Fig. [5l results of S(k) at the same densities as in 
Fig. [H are reported. The results again show the ex- 
pected behavior: with the increase of p, the strength of 
the main peak increases and moves to higher momenta in 
a monotonic way. At low momenta, the slope of S(k) de- 
creases with the density, following the limiting behavior 
linifc^o S(k) = hk/(2mc) driven by the speed of sound c. 
The DMC data start at a finite k value inversely propor- 



tional to the box size L. 

A characteristic signature of bulk superfluidity is the 
finite value of the condensate fraction uq, i.e. the fraction 
of particles occupying the zero-momentum state. It has 
been extracted from the long-range behavior of the one- 
body density matrix p(r), by means of the asymptotic 
condition tiq = lim r ^oop(r). We have verified, by in- 
creasing the number of particles of the simulation at dif- 
ferent densities, that the size dependence of no is smaller 
than its statistical error. The results obtained using the 
extrapolated estimator uq ~ 2(no)DA/c ~ ( n o)vMC are 
presented in Fig. [51 The line on top of the data corre- 
sponds to the exponential fit 

no(p) = Aexp(-bp) (9) 

with A=3.6(3) and 6=455(11) A -1 , which describes well 
the DMC data. At the equilibrium density, no=0. 129(3). 
For comparison, 4 He at the equilibrium has no=0. 084(1), 
as obtained by the DMC method using a Jastrow wave- 
functions for importance sampling in Ref. [HI 

Previous theoretical estimates of the energy per par- 
ticle in bulk tritium were obtained by VMGii or 
Brueckner-Bethe-Goldstone formalism. 12 The TJ-TJ in- 
teraction was modeled with the Morse potential and the 
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calculations were performed for densities up to 0.015 
A~ 3 . Etters et al^ obtained an equilibrium energy per 
particle E/N — —0.75(2) K for the equilibrium density 
po = 0.0463 A~ 3 . For the same density, Joudeh et al^ 
obtained —0.759 K. Comparing this equilibrium density 
with the present DMC results one notices that this vari- 
ational estimate lies below our estimated spinodal point. 
In order to compare our results with previous ones, we 
have carried out simulations at a density p — 0.005 A -3 
using both the Morse potential and the JDW one. With 
the Morse potential and at the VMC level we obtain 
E/N = —0.74(1) K, which is within the errorbars the 
same result as that of Etters et al.(— 0.75(2) K). Using 
the Morse potential, the DMC calculation lowers the en- 
ergy to -0.978(1) K. At the same density 0.005 A~ 3 , 
we have also used the JDW potential obtaining at the 
DMC level a sizeable lower energy of -2.900(2) K. The 
observed differences in the equation of state imply also 
substantial changes in the estimated pressure and com- 
pressibility. For example, for a density p = 0.01 A -3 
we find a pressure of 11.6(1) bar, while in Ref. [ll| it is 
estimated to be 17 bar. At higher densities, the differ- 
ence between the pressure results grows, amounting to 
approximately 80 bar at 0.015 A -3 . 

Using the relation for the transition temperature of the 
ideal Bose-Einstein gas 

r -= 3 - 3i (i£K' (io » 

Etters et ali^- estimated the temperature of superfluid 
transition of spin-polarized T to be 1.48 K. With the 
same argument, our results for the equilibrium density 
suggest a superfluid transition temperature of 2.02 K. 
This is of course only an approximate estimation since in 
liquid 4 He Eq. fTTJl gives 3.1 K instead of the right one of 
2.17 K. 



B. Solid phase 

We have performed calculations of the spin-polarized 
T solid phase with three different lattices (bcc, fee, hep) 
and using the Nosanow- Jastrow wave function §5§ for im- 
portance sampling. The energy per particle in bcc, fee, 
and hep lattices has been obtained for different densities 
in the range from 0.008 A~ 3 to 0.024 A" 3 . Our DMC re- 
sults show that the energies per particle for the three lat- 
tices are statistically indistinguishable in all the studied 
density regime. Still, it reaches the lowest values in the 
hep solid phase, so this lattice seems to be energetically 
preferred. As an example, the results at two densities for 
bcc, fee and hep lattices are respectively: for p = 0.011 
A~ 3 , near the melting density, the energies per particle 
are E/N = -1.96(5), -1.98(7), and -2.04(6) K; at a 
higher density p = 0.018 A" 3 , E/N = 15.8(2), 15.3(2), 
and 15.26(8) K. The same behavior has been observed for 
all densities and therefore we decided to investigate solid 



P (A" 3 ) 


E/N (K) 


T/N (K) 


P (bar) 


c (m/s) 


0.01 


-2.56(6) 


14.16(8) 


4.7(7) 


321(18) 


0.011 


-2.04(6) 


16.51(9) 


11(1) 


404(19) 


0.015 


4.17(8) 


27.89(11) 


83(4) 


782(28) 


0.018 


15.26(8) 


38.30(12) 


218(8) 


1113(36) 


0.024 


62.59(8) 


62.54(13) 


902(29) 


1898(58) 



TABLE II: Results for solid Tl at different densities p: en- 
ergy per particle (E/N), kinetic energy per particle (T/N), 
pressure (P), and speed of sound (c). Figures in parenthesis 
are the statistical errors. 

TJ, properties assuming its hep crystalline structure. It 
is important to notice that the bcc lattice proved to be 
energetically preferred in a recent study of the gas-solid 
phase transition in HJ.— as well as in a study of solid 
hydrogen at very high pressure^. 

In Fig. [3 the DMC energies per particle of the solid 
phase at different densities have been shown, for the three 
lattices. The line on top of the data in the figure corre- 
sponds to the equation of state of the solid hep lattice. 
It has been obtained by fitting the DMC results with a 
polynomial function of the form 

e(p) = s 2 p 2 + s 3 p 3 + s 4/ o 4 , (11) 

with parameters s 2 = -10.47(11) x 10 4 KA 2 , s 3 = 
7.13(15) x 10 6 KA 3 , and s 4 = 7.3(5) x 10 7 KA 4 . 

The total and kinetic energies per particle for several 
selected densities are given in Table |TTJ Kinetic energy 
is determinated in the same way as in the liquid phase. 
The total energy is neg ative for p < 0.012 A" 3 , while 
for greater densities the kinetic energy exceeds the abso- 
lute value of the potential one causing the total energy 
to become positive. The potential energy is negative for 
all the considered densities, with a single exception cor- 
responding to the highest density 0.024 A~ 3 . As it is 
reported in Ref. [l|, the smaller mass of HJ, atoms caused 
the potential energy to enter in regime of positive values 
for slightly smaller densities (p > 0.02 A -3 ). 

In FigJSJ we show the pressure and speed of sound of 
solid TJ, obtained from the equation of state Ijlip using 
the thermodynamic relations (0 and §8§ as a function of 
the density. Comparison with values of the same quanti- 
ties in solid bulk HJ,, reported in Ref. Q], reveals smaller- 
pressure and speed of sound for solid TJ, at the same 
densities. 

The spatial pattern of the solid structure is reflected 
in the two-body radial distribution function g(r). In Fig. 
U3 we have plotted g(r) for same selected densities. The 
strength of the main peaks are greater in the solid than 
in the liquid phase, as can be seen by comparison of our 
results at densities 0.01 A~ 3 and 0.02 A -3 . Also, it is 
clear that secondary peaks are larger in the solid phase. 
Just like in the liquid, when the density increases the 
height of the main peaks grows and moves to shorter 
distances. 
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FIG. 7: Energy per particle of solid TJ. as a function of the 
density p for the hep (boxes), fcc(triangles) and bcc(circles) 
lattices. The solid line corresponds to the fit to the DMC 
energies using Eq. . The error bars of the DMC energies 
are smaller than the size of the symbols. 
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FIG. 9: Two-body radial distribution functions of the solid 
phase. From bottom to top in the height of the main peak, 
the results correspond to densities 0.01, 0.011, 0.015, 0.02 and 
0.024 A" 3 . 
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FIG. 8: Pressure and speed of sound of solid TJ. as a func- 
tion of the density. Left (right) scale corresponds to pressure 
(speed of sound). 
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FIG. 10: Maxwell construction based on plotting the energy 
per particle, E/N as a function of 1/p. The densities at which 
the first-order transition occurs are identified by finding the 
common tangent (solid line) to both the solid (dotted line) 
and liquid curve (dot-dashed line). 



C. Liquid-solid phase transition 

An important information that can be derived from the 
DMC simulations of the liquid and solid phases is the lo- 
cation of the liquid-solid phase transition point. As in 
a recent investigation of the gas-solid phase transition 
in H|— , we have used the double-tangent Maxwell con- 
struction to determine the transition. This well-known 
method is based on the search of a common tangent to 
both the liquid and solid equations of state whose inter- 
sections give the freezing (pf) and melting (p m ) densities, 
as plotted in Fig. [TDJ Using the equation of state of the 
hep crystalline structure, we have obtained pf = 0.00964 
A" 3 = 0.477 a" 3 and p m = 0.01069 A~ 3 = 0.528 a' 3 , 
which corresponds to a common pressure at the phase 
transition of P = 9(1) bar. In addition, in order to 
estimate the influence of the lattice type on the tran- 
sition pressure, we have also calculated the transition 
densities and pressure for fee and bee lattices. For a 
fee lattice P = 9.5(1.0) bar, while for bec we have ob- 
tained P = 9.9(1.0) bar. The three transition pres- 
sures fall within the errorbars, but since hep is consis- 
tently lower than the others it leads us to conclude that 



it may be slightly preferred, like it is preferred in 4 He. 
The transition densities in 4 He are pi = 0.430 <r~ 3 and 
p m = 0.468 a~ 3 ^2. corresponding to a pressure of 25.3 
bar. The comparison can also be made with with gas- 
solid transition in spin-polarized HJ, 1 . In this case, the 
phase transition occurs at higher densities pf = 0.01328 
A~ 3 , p m = 0.01379 A~ 3 and much higher pressure, 
P = 173(15) bar, than in spin-polarized T. This effect 
can be explained as a consequence of its isotopic dif- 
ference since TJ. atoms have approximately three times 
greater mass than HJ. atoms. 

As in any first-order phase transition, the density is 
discontinuous in the transition point. Another quantity 
which is also not continuous crossing the phase transi- 
tion is the kinetic energy per particle. Namely, near the 
freezing density the kinetic energy per particle of the liq- 
uid is T/N = 13.06(3) K, while at the same time, near 
the melting density but in the solid phase the same en- 
ergy is T/N = 16.51(3) K. Therefore, our results for bulk 
TJ. show a discontinuity of the kinetic energy of around 
3.5 K; in bulk HJ., the same difference has been shown to 
be more than twice times larger ii On the other hand, the 
value of the condensate fraction in the two spin-polarized 
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FIG. 11: Two-body radial distribution function at the liquid- 
solid phase transition. The solid line corresponds to the liquid 
at pi and the dashed line to the solid at p m . 



systems at the corresponding transition densities are sim- 
ilar, 0.03 for TJ. and 0.04 for H|-i. 

A measure of the mean displacement of TJ, atoms with 
respect to the lattice sites is obtained by computing the 
Lindemann's ratio (7). As usual, we calculate it by sam- 
pling the expression 7 = yj ((r — ri) 2 )/a,L, where is 
the lattice constant. The parameter 7 decreases mono- 
tonically with increasing density and the values we have 
obtained show a slight dependence on the particular lat- 
tice chosen for the simulation: it is greatest for hep 
and smallest for bcc. Particularly, at the melting den- 
sity of the hep phase Lindemann's ratio assumes a value 
7 = 0.26, which is the same as in He and similar to the 
value estimated for solid HJ (7 = 0.25)i. 

The discontinuity in the liquid-solid transition is also 
revealed in the difference between g(r) of the liquid and 
solid phases (Fig. Hip . Even more dramatic difference 
between the two phases is demonstrated in Fig. 1 121 where 
DMC results of S(k) are shown at pf and p m . S(k) 
in the solid phase is characterized with strong peaks at 
reciprocal lattice sites whereas this behavior is clearly 
not observed in the S(k) of the liquid phase. Finally, 
it is worth noticing that the main peak of solid S(k)is 
slightly weaker in TJ than in HJ (reported in Ref. [0). 
The main reason for this lies in the fact that the liquid- 
solid transition in TJ emerges at lower densities than 
gas-solid transition in HJ,. 

IV. CONCLUSIONS 

The ground-state properties of spin-polarized tritium 
TJ, have been accurately determined using the DMC 
method in both the liquid and solid phases. The ob- 
tained results are based on the precise knowledge of the 
TJ.-TJ, interatomic potential, which combined with the 
accuracy of the DMC method, allowed for a nearly ex- 



act determination of the main properties of the system. 
All the previous results on liquid TJ, were based on a 
Morse potential and approximate computational meth- 
ods. Hence, our predictions for equilibrium density and 
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FIG. 12: Static structure factor at the liquid-solid phase tran- 
sition. The results correspond to the liquid at p{ and to the 
solid at pm. 

energy per particle differ substantially from those. We 
find that the equilibrium density is p = 0.00747(1) A~ 3 , 
which expressed in units of <r~ 3 is very similar to the one 
of liquid 4 He (p = 0.369 <r~ 3 for TJ^ p = 0.365 cr" 3 in 
4 He ). Previous predictions of the equilibrium density lie 
below our predicted spinodal density. Despite a similar 
equilibrium density, TJ has approximately half the equi- 
librium energy per particle of liquid 4 He as a consequence 
of its smaller mass and shallower potential. 

At a high enough density the system freezes. We 
have studied the energetic and structural properties of 
the solid phase for three lattices. Differences in ener- 
gies between the phases are almost indistinguishable, but 
hep seems to be slightly preferred over the fee and bcc 
ones. From the DMC equations of state of the liquid and 
solid phases, we have localized the liquid-solid transition 
point of TJ, for the first time, to the best of our knowl- 
edge. At zero temperature, the phase transition occurs 
at P = 9(1) bar. 
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